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We analyze a system of globally coupled logistic maps with asynchronous updating. We show that 
its dynamics differs considerably from that of the synchronous case. For growing values of the 
coupling intensity, an inverse bifurcation cascade replaces the structure of clusters and ordering in 
the phase diagram. We present numerical simulations and an analytical description based on an 
effective single-element dynamics affected by internal fluctuations. Both of them show how global 
coupling is able to suppress the complexity of the single-element evolution. We find that, in contrast 
to systems with synchronous update, internal fluctuations satisfy the law of large numbers. 



PACS: 05.45-fb, 05.90-fm 
I. INTRODUCTION 

Introduced by Kaneko in 1984 coupled map lat- 
tices have proved to be a powerful tool in the study of 
spatioternporal chaos and pattern formation in complex 
systems l^-Ql- They have found applications as mod- 
els in different branches of science Globally cou- 
pled maps 1^, which bear similarity to the Sherrington- 
Kirkpatrick model, constitute a kind of mean-field exten- 
sion of coupled map lattices. In an ensemble of globally 
coupled maps, the elements are seen to form clusters of 
synchronized activity. As the coupling increases, a vari- 
ety of phases can be identified — from coherent, through 
ordered and partially ordered, to turbulent — depending 
on the number of clusters in the system. For sufficiently 
large coupling intensities, the generic behavior consists 
of full synchronization of the whole population, with all 
the elements having identical evolution. 

In usual models of globally coupled maps the elements 
update their state at the same time, namely, their evo- 
lution is synchronous. For an ensemble of N elements 
whose individual dynamics is given by the nonlinear map 
Xi{t + 1) = f[xi{t)], a typical scheme for global coupHng 
is given by 



N 



x,{t + !) = (!- e)f[xdt)] + /[^.- W]' * = 1, • • • , ^, 



(1) 



where e £ [0, 1] is the coupling intensity. Equation 
is applied synchronously to all the maps of the system, 
with the values of all the sites at the previous time as 
inputs. 

From a realistic viewpoint, however, synchronous up- 
dating does not seem to be very plausible in models of 
real systems 0-|l5|. As discussed, for instance, for neu- 
ral networks , an independent choice of the times at 



which the elements of a given complex system update 
their respective states should provide a closer approxi- 
mation to reality. It has however to be mentioned that, 
in systems with local coupling, asynchronous update has 
been shown to lead to trivial behavior We present 
here a first analysis of asynchronous globally coupled 
maps and find their behavior to be completely different 
from that of the usual synchronous models. We have 
analyzed two possible asynchronous schemes: 

(a) At each time step, update the elements according to 

the sequence given by a prescribed "list" that con- 
tains every element once. The order of the elements 
in the list is chosen at random at each step. 

(b) Divide each time step into N substeps. At each sub- 

step, update one element chosen at random. After 
N substeps the map operates, on average, once over 
each element in the system. 

For the sake of clarity, let us formulate mathematically 
the second scheme. Let ri{t) be a stochastic process that 
takes an integer value ry — 1,2, ... ,N at each substep, 
with uniform probability. The dynamics proceeds in the 
following way: 



il~e)f[x,{t)] 

Xi{t) 



e 

N 



if * = 77(0, 
(2) 



In the first scheme, rj{t) takes the values 1,2,...,N at 
random, but only once within each time step. It can 
be seen that, in this dynamics, the mean field F{t) — 
N~^^j f[xj{t)] has different values when acting over 
different elements, as time proceeds by steps of length 
This feature and the fact that asynchronous updat- 
ing incorporates a stochastic component in the evolution, 
contrast with the deterministic synchronous scheme M. 
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Evolution under the two asynchronous schemes intro- 
duced above displays the same global properties. How- 
ever, in the second scheme one or more elements may 
fail to be updated in some time step. As a consequence, 
since N updates have to be made during a time step, 
some other elements will be updated more than once. 
This has nontrivial consequences in the dynamics, to be 
discussed below. Besides, the second scheme has the ad- 
vantage of being much faster computationally. 

In this paper, we analyze an ensemble of globally cou- 
pled logistic maps. We show that the dynamics under 
asynchronous updating substantially differs from that of 
the synchronous case. The structure of clustered phases 
completely dissapears. In its place, an inverse bifurcation 
cascade develops as the coupling intensity is increased. 
In the following, we analyze this phenomenon from nu- 
merical simulations and construct a phase diagram for 
the different regimes displayed by the dynamics. Then, 
we propose an analytical explanation for the bifurcation 
cascade, whose results compare successfully with the nu- 
merical data. Finally, we discuss how the fluctuations 
generated in the internal dynamics are reflected in the 
average evolution. 
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FIG. 1. Bifurcation diagrams of a system of A'^ = 1000 
logistic maps, as a function of the coupling intensity e. The 
diagrams are constructed following the evolution of a single 
element taken at random, with (a) A = 3.5, (b) A = 3.84, (c) 
A = 4. 



II. INVERSE BIFURCATION CASCADE AND 
PHASE DIAGRAM 

In this section we present extensive numerical simula- 
tions of the evolution given in Eq. for the standard 
logistic map, f{x) = Xx{l — x). We concentrate in the 
values of A where this map displays its bifurcation cas- 
cade, 3 < A < 4, leading through period doubling from a 
stable fixed-point state to a completely developed chaotic 
regime. For reasons that will become evident later, we re- 
strict the coupling intensity e to the interval < e < 0.5. 

In Fig. |l| we show the bifurcation diagrams for the 
evolution of a single element, chosen at random from a 
population of 10"^ maps, as a function of e, for three val- 
ues of A. These diagrams are constructed as usual. For 
given values of e and A the system is let to evolve until 
transients have elapsed. Then, the state of the chosen 
element is recorded and plotted at some successive time 
steps. 

Figure |l|(a) shows the bifurcation diagram for A = 3.5, 
where the (uncoupled) logistic map evolves in a period-4 
orbit. As e grows, the four initial branches collapse into 
two branches which, in turn, merge into a single branch. 
Increasing the coupling intensity leads thus the evolution 
of a single element to display an inverse bifurcation cas- 
cade. The same scenario occurs for other values of A. 
For A = 0.38, where the logistic map is within its largest 
period-3 stability window, increasing e leads the evolu- 
tion to successively exhibit chaotic regimes of one and 
two bands and, eventually, stable branches that finally 
collapse into a single stable state (Fig. |l|(b)). For A = 4, 
in the extreme chaotic regime of the logistic map, the full 
bifurcation diagram is run over backwards (Fig. |l|(c)). 

It is apparent from Fig. ^ that the evolution of a sin- 
gle element in the ensemble is subject to the action of 
noise. Indeed, the branches of nonchaotic evolution are 
not perfectly defined — except in the case of a single sta- 
ble state — and high-order bifurcations, close to the onset 
of chaos, are clearly suppressed. This noise, which is 
produced internally in the system via the asynchronous 
updating, induces spreading of the states visited during 
the evolution when the coupling is different from zero. 

The branches of nonchaotic evolution, with their small 
but noticeable dispersion, are what remains of the stable 
periodic orbits of the deterministic map under the effect 
of the internal noise. In these branches, the evolution 
driven by the two schemes presented in the Introduc- 
tion differ. Scheme (a) defines true noisy periodic orbits, 
where each element visits the set of available states fol- 
lowing the same periodic sequence as the deterministic 
map. On the contrary, scheme (b) allows some elements 
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to skip the update in a time step. Consequently, some 
other elements are updated more than once. As a result, 
the observed successive states do not follow a periodic 
sequence. To simplify the discussion, in the following wc 
will call these regimes of noisy nonchaotic (periodic or 
nonperiodic) evolution "period-rt" orbits, referring rather 
to the set of states visited by each element. 

The chaotic bands, meanwhile, are also blurred by the 
internal noise. We will thus refer to a regime where a 
noise-free map would exhibit deterministic chaotic evo- 
lution as an "n-band" chaotic region. 

To summarize the variety of behavior observed for a 
single element as e and A vary, we have constructed a 
phase diagram in the plane spanned by these two param- 
eters for a population of 10'^ elements. For each value of 
e and A the evolution has been calculated with scheme 
(b) during 10* time steps. After a transient of 5 x 10'^ 
steps, a 500-column histogram over the states visited by 
a single element in the remaining 5 x 10"^ steps has been 
produced to identify the kind of evolution corresponding 
to those values of the parameters. The phase diagram is 
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FIG. 2. Phase diagram of a system of A*' = 1000 logistic 
maps. The two lines are the analytic curves given by Eqs. (^) 
and (|). 

The large upper-left region corresponds to a fixed 
point. In this zone of the parameter space, from any ini- 
tial condition, all the elements are attracted to the same 
stable fixed point. Without coupling — i.e. for e = 0, 
on the horizontal axis of the plot — such behavior would 
be observed for A < 3 only. At larger values of A and, 
notably, even up to the fully developed chaotic regime 
(A = 4), there always exists a coupling intensity able to 
suppress the complex behavior of the uncoupled system. 
As mentioned above, in this fixed-point regime internal 
noise is also suppressed. 



Below the fixed-point region, there is a series of zones 
shown, alternately, with black squares and empty space 
in Fig. 1^. In this phases the evolution of single elements 
display "periodic" and "chaotic" noisy orbits, as indi- 
cated in the plot. Due to the noise, the limits between 
these regions are not sharply defined and some of the 
zones — such as the period-8 and the period-5 regions — 
are truncated at certain values of A and e. Between the 
zones of period-4 orbits and two-band chaotic evolution, 
higher-order periodic orbits as well as chaotic evolution 
in more than two bands are missing. The period-3 stabil- 
ity window is instead clearly detected, immersed in the 
one-band chaotic regime. 

Since blurring of boundaries between zones and sup- 
pression of higer-order periodic and chaotic orbits are a 
direct consequence of the internal noise, it is expected 
that — if noise decreases upon incresing the number of 
elements in the ensemble — the phase diagram becomes 
more sharply defined for higher N. It is however known 
that, in synchronous globally coupled maps, the inter- 
nal noise does not decrease in the usual self-averaging 
way as N~'^/'^ — as predicted by the so-called law of large 
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numbers — but in a much slower manner |17 
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tion IV we analyze this aspect for asynchronous evolu- 
tion, finding that the law of large numbers is recovered 
in these systems. 

Note finally that if, in contrast with Fig. |l|, the bifurca- 
tion diagram would be plotted for fixed e as a function of 
A, the bifurcation cascade would proceed forwards. For 
large values of the coupling intensity, however, higher- 
order bifurcations and chaos would not be reached even 
for A = 4. 



III. MEAN FIELD APPROACH 

It is possible to describe analytically some features of 
the bifurcation cascade by performing a kind of mean 
field approximation to Eq. (||). From the viewpoint of 
a single element, updating of its state occurs, in aver- 
age, once per time step. At the successive times where a 
given element is updated, the sequence of values of the 
mean field F{t) = N-^ J2j fi^jit)] fluctuates due to the 
evolution of the individual states of all the elements, ef- 
fectively resembling a stochastic process. In such a way, 
the single-element dynamics can be though of as given by 
a deterministic map subject to the action of an effective 
"external" stochastic forcing. 

In order to characterize this effective forcing, we 
assume — as suggested by the numerical simulations — 
that the evolution of the system (||) determines, at long 
times, a well-defined measure on the space x. At 
any time step the value of the elements will distribute 
according to this measure and, for N oo, the mean 
field F{t) will approach a constant Fq = J f{x)ii{x)dx. 
For finite N , however, the field F{t) will fluctuate around 
Fq, so that we can write 
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lim F{t)^Fo+m, 



(3) 



where ^{t) is a stochastic process of zero mean. 

Within this approximation, the dynamics of a single 
element is 



x,{t + !) = (!- e)f[x^{t)] + eFo + e^t), 



(4) 



The effective individual evolution is thus given by a de- 
terministic renormalized map = (1 — e)f{x) + ei^o 
upon which an additive noisy force of zero mean acts. It 
has to be noticed that the constant Jq, being determined 
by the collective dynamics, prevents the population from 
complete decoupling. 

It is well known that a deterministic map subject to 
the action of a moderate additive noise of zero mean pre- 
serves its critical behavior, the only effect of noise being 
the suppression of high-order bifurcations. Within this 
picture, we can then cxplicitcly calculate the first bifur- 
cation points as follows. 

Suppose firstly that the map /"^{x) has a stable 
fixed point xq as its only attractor. Since — in the ab- 
sence of noise — all the elements are attracted to it, we 
have fj,{x) = 6{x — xq) and Fq = f{xo)- The effec- 
tive single-element evolution becomes thus Xi{t + I) = 
(1 — €)f[xi{t)] + ef{xo). The equation 



2^0 = (1 - ()f{xo) + ef{xQ) = /(a;o), 



(5) 



determines xq, which results to coincide with the fixed 
point of the original map f{x). The evolution equation 
can be now readily linearized to find that xq is stable if 
the coupling intensity satisfies 



i-i/'(^o)r 



(6) 



For the logistic map, that has a stable fixed point at 
xo = 1 - A^^ for 1 < A < 3, Eq. (||) implies that the 
state Xi — Xq is stable if e > 1 — |2 — A|~^. That is, for 
2 < A < 3 the state is stable for all e and, for A > 3 it is 
stable for e large enough. In Fig. ^ we have drawn the 
curve el{X) that defines the boundary between the fixed 
point and the period-2 regions. This analytical result is 
in full agreement with simulations. 

A similar argument can be applied to period-2 orbits. 
Due to the randomness of the updating scheme, at each 
time step practically one half of the elements will be 
at one of the two states of the orbit, say xa, and the 
other half at the other, xb- The corresponding invari- 
ant measure is /Lt(a;) = [S{x — xa) + S{x — xb)]/'^ and 
^0 = [fi^A) + /(a^B)]/2. Therefore, the equations which 
determine the values of xa and xb are 



XA - {l-e)f{xB) + e[f{xA) + f{xB)]/2 
XB = {l-e)f{xA)+e[f{xA)+f{xB)]/2. 



(7) 



These equations are equivalent to xa = /{xb) and 
Xb = /(xa), so that Fq = (xa + xb)/2 and xa,b are 
also the two states corresponding to the period-2 orbit of 



the original map f{x). For the logistic map these equa- 
tions can be solved, and the stability condition implies 



e > e,. = 



-1 - 2A + A^ - VI -6A + 3A2 
A(A-2) ■ 



(8) 



The line e^W is also shown in Fig. H, separating the re- 
gions of period-2 and period-4 orbits. This result is again 
in good agreement with simulations, though a slight de- 
viation is observed for large values of e, i.e. near A = 4. 
This deviation can be ascribed to the increasing effect of 
noise for growing e (cf. Eq. (^)), that blurs the bound- 
aries where bifurcations occur. 

Analogous reasoning may be used to determine the 
boundaries between other zones. For higher-order pe- 
riodic orbits, however, the equations cannot be explicitly 
solved in the case of the logistic map. For chaotic evo- 
lution, moreover, the measure fi{x) should be obtained 
numerically. 



IV. FLUCTUATIONS OF THE MEAN FIELD 

A key ingredient in the analysis presented in the pre- 
vious section is the assumption that, for large popula- 
tions, the system (|^) defines a measure in the one-element 
state space x, such that for sufhciently long times the 
elements exhibit a well-defined distribution in x. It is 
implicit in that assumption that the amplitude of inter- 
nal noise decreases as A'' grows. As mentioned above, 
however, previous work has shown that fluctuations in 
deterministic synchronous globally coupled systems does 
not self-average in the usual way ||l^,|l^. It is therefore 
worthwhile to analyze this point in some detail for the 
present case of asynchronous update. 

The fluctuating part of the mean field, ^(t) in Eq. 
(^), can also be analytically studied within some fur- 
ther approximations. In order to illustrate our argu- 
ments, let us restrict ourselves to the period-2 orbit, 
where fi{x) = [S{x — xa) + S{x — xb)]/^. 

Due to the fluctuations arising at finite values of N, 
the number of elements in each state will not be A^/2, but 
rather a fluctuating number n^(i) at xa and nsit) = N— 
nA{t) a,t Xb- Since the mean field F{t) will consequently 
differ from Fq = [/(xa) + fixB)]/^ = {xa + xb)/2, this 
implies that the individual states will spread around the 
values XA and xb- At a given time, therefore, the ele- 
ments near xa.b will have states ^ = XA.B+Sxi^ ^ (t). 
The fluctuating mean field will thus be 



F{t)^ 



1 

N 



^ f{xA + Sxi^) + ^ f{xB + Sxis) 



(9) 



where the sums Vl. and Y^,- run over the elements near 
XA and Xb, respectively. 

In the lowest-order approximation, we have F{t) = 
[uAfixA) + nBf{xB)]/N = {uaxb + nBXA)/N that can 
be rewritten as 
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Fit) 



XA + Xb 



N 



nA{t) 



N 



{xb~xa). (10) 



In this expression we identify the fluctuations of F{t) 
as ^{t) = [nA{t) - N/2]{xb - Xa)/N. Since updat- 
ing is apphed at random in the ensemble, nA{t) is ex- 
pected to have a binomial distribution around iV/2. This 
readily shows that £^{t) is a stochastic process with zero 
mean, whose mean square dispersion depends on N as 
^ 7V~^/^. We also note that the amplitude of fluctua- 
tions depends Hnearly on the difference jxyi — xb] between 
the two states of the period-2 orbit. As a byproduct, this 
indicates that the fixed-point state exhibits no fluctua- 
tions. In Fig. H we have plotted A^ as a function of 
\xa — xb\ for a fixed value of A in a population of 10"* 
logistic maps. The coupling intensity e runs over the re- 
gion of period-2 orbits, making | — a^s | vary from small 
values near t = e\io larger values as e decreases. A large 
region in which A^ varies linearly with — a;_B| is clearly 
observed. 
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FIG. 3. The standard deviation of the fluctuations of the 
mean field, A^, as a function of the the separation between the 
two states of the period-2 orbit. The system size is N = 10* 
and A = 3.4; e approaches from below. The abcisas used to 
plot the data correspond to the values of xa and xb provided 
by Eq. (Q) as a function of e. 

Remarkably, the proportionality of the mean field fluc- 
tuations with iV^^/^ is also numerically verified in the re- 
gions with more complex dynamics. The validity of the 
law of large numbers proves thus to be a generic feature in 
the fluctuations of the asynchronous ensemble of globally 
coupled maps. In Fig. ^ we show a double logarithmic 
plot of the mean square deviation AF of the mean field as 
a function of the system size for e = 0.1 and A = 4, i.e. in 
the chaotic regime. The A^~^/^-dependence is apparent. 
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FIG. 4. The standard deviation of the mean field, AF, as 
a function of the system size. The double logarithmic plot 
shows the decaying behavior N~^^^ for system sizes ranging 
from TV = 300 to 10*^ for A = 4 and e = 0.1. A hne of slope 
— 1/2 is also shown for comparison. 



V. MEAN FIELD DEPENDENCE ON THE 
COUPLING INTENSITY 

In order to complete our discussion, it is worthwhile 
to investigate how the collective behavior varies as the 
coupling intensity is modified. In this section we present 
numerical results on the dependence of F{t) on e for a 
given value of A in an ensemble of logistic maps. 

In Figure ^ we show the time-averaged mean field (F) 
as a function of the coupling intensity, for A = 4 and 

= 10"^. At e = the elements are completely uncou- 
pled. In the fully chaotic regime (A = 4), the individual 
states are symmetrically distributed in the interval (0,1) 
according to the invariant measure of the logistic map. 
This uncorrelated state is then characterized by an aver- 
age field (F) = 0.5. For nonzero values of the coupling 
intensity e, we observe that a correlated state develops, 
characterized by values of (F) larger than 0.5. The aver- 
age mean field grows however in a highly nonmonotonic 
way. Large downward peaks can be observed, where the 
collective state again approaches the value (F) — 0.5. 
The broadest of these peaks, near e — 0.07, coincides 
with the period-3 stability window. The inset of Fig. |^ 
shows a detail of the same function in the region of small 
coupling intensity (for an ensemble of 5 x 10^ elements). 
The average mean field exhibits a complex dependence 
on e, seemingly, with self-similar features. These features 
are probably reflecting the presence of higher-order sta- 
bility windows in the chaotic regime. For large values of 
e the system abandons the chaotic state and, for e « 0.2, 
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it enters the regime of periodic orbits. Correspondingly, 
(F) behaves less erratically. 



Abramson and Zanetle. Figure 6. 
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FIG. 5. Temporal average of _F as a function of e. The 
system has A'^ = 1000 elements and A = 4. Inset: a detailed 
view of the region of small e, with A'^ — 50000. 

For lower values of A a similar picture can be observed. 
In such cases, of course, the average mean field does not 
start at = 0.5 for e = 0, since the values of Xi do not 
cover the interval (0, 1) with a symmetric distribution. 
As the coupling intensity grows, at the critical value ej, 
the entire population collapses into the fixed point and 
(F) reaches a fixed maximal value. 

Finally, we have studied how the fluctuations of F{t) 
around its average value (F) depend on the coupling in- 
tensity. In Figure |^ we have plotted the mean square dis- 
persion of F, A_F = ^ {(F — (F))^), as a function of e, for 
several system sizes. These plots show a rather uniform 
background, that disappears at e = 0.5 when the whole 
system is attracted to the fixed point (the non-linear pa- 
rameter is A = 4). This fiuctuation level is reduced by 
enlarging the system. Superimposed to this background, 
some sharp spikes of enhanced fluctuations are seen for 
low values of e. They coincide with the downward peaks 
of (F) in Figure ||, and correspond to the stability win- 
dows in the chaotic regime. As illustrated in the inset 
of Fig. ^ for the widest period-3 window, the behavior 
of any element of the ensemble in these regions is highly 
intermittent. During certain time intervals, the elements 
are engaged in periodic orbits but, occasionally, they ex- 
hibit a regime of chaotic evolution. This intermittency 
between two qualitatively different forms of motion, each 
of them having specific values of (F) and AF, causes 
the overall mean-field dispersion to attain unusually large 
levels. These anomalous fiuctuations can even grow upon 
increasing the system size. 
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FIG. 6. The mean square dispersion AF as a function of 
e, for several system sizes. Dotted line: N ^ 10^; full thick 
line: iV = 10*; full thin line: N ^ 10^. Increasing N, the 
background is seen to diminish, while the peaks are seen to 
grow. For N — 1000, the peaks cannot be resolved from the 
background. Inset: The orbit of a single element at A = 4 
and e — 0.069, from a system with 1000 elements, displaying 
intermittent behavior between a period-3 orbit and chaotic 
motion. 



VI. SUMMARY AND CONCLUSIONS 

In this paper we have analyzed the collective behav- 
ior of an ensemble of globally coupled maps whose dy- 
namics is updated asynchronously. Asynchronous up- 
dating, whose study is motivated by the aim at realisti- 
cally modelling the evolution of real systems, introduces 
a stochastic ingredient in the dynamics, with nontriv- 
ial consequences in the behavior of the population. Al- 
though some of our analytical results hold for any kind of 
coupled maps, we have focused our attention — in partic- 
ular, in the numerical simulations — in the case of logistic 
maps, f{x) = Ax(l — x). 

We have numerically found that, for a fixed value of the 
parameter A, increasing the coupling intensity leads the 
system to simpler and simpler evolution, running back- 
wards over the bifurcation diagram of the logistic map. 
For sufficiently large coupling, in fact, a fixed-point state 
is reached for the whole ensemble, even when A corre- 
sponds to chaotic individual motion. This is in strong 
contrast with the behavior of globally coupled maps with 
synchronous updating. Indeed, large coupling intensities 
lead such systems to a synchronized state where all the 
elements reproduce the evolution of a single, uncoupled 
map. Below the synchronization threshold, those sys- 
tems exhibit a regime of clustering not observed in the 
present case of asynchronous updating. 

An approximate analytical description for the asyn- 
chronous ensemble can be achieved by constructing an ef- 
fective single-element dynamics, which result to be driven 
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by the map f {x) = (1 — e)/(a;) + eFo, where e is the cou- 
phng intensity and Fq is a constant to be determined self- 
consistently from the collective evolution. For the case of 
logistic maps, this approximate picture makes it possible 
to analytically calculate the threshold of the lower-order 
bifurcations. These results compare successfully with nu- 
merical simulations. Note that as e increases, the weight 
of the nonlinear term in the effective dynamics decreases, 
explaining why the bifurcation diagram develops back- 
wards when the coupling intensity grows. 

The effective individual dynamics is affected by inter- 
nal fluctuations, which enter the single-element evolution 
as an additive noise term. As is well known for noisy 
maps, the main effect of these fluctuations consists of 
suppression of higher-order bifurcation and blurring of 
both regular and chaotic motion. We have studied how 
the amplitude of the internal noise depends on the en- 
semble size, i.e. on the number N of elements in the 
population, by analyzing the temporal behavior of the 
mean field F{t) = J2i f[^i{t)]- Numerical simula- 
tions show that the fluctuations of F{t) decreases with 
the system size as N~^/^. This is again in strong contrast 
with synchronous globally coupled maps, for which it is 
known that fluctuations decrease in a much slower man- 
ner. Here, instead, the intrinsic stochastic character of 
the evolution makes fluctuations to obey the usual selfav- 
erage statistics, and the law of large numbers holds. This 
result suggest that the violation of the law of large num- 
bers is not a robust feature upon introduction of random 
elements in the dynamics of globally coupled ensembles. 

In summary, we have shown that the collective behav- 
ior of globally coupled maps with asynchronous updating 
exhibits important differences when compared with that 
of synchronous dynamics. Coupling is able to suppress 
the complexity of individual evolution and internal fluc- 
tuations selfaverage in the usual way as the system size 
increases. The extension of these results to ensembles 
formed by more complex maps and by continuous-time 
dynamical systems should be the subject of further anal- 
ysis. 
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